Comparative population genomics in Tabebuia alliance shows evidence of adaptation in Neotropical tree species.

The role of natural selection in shaping spatial patterns of genetic diversity in the Neotropics is still poorly understood. Here, we perform a genome scan with 24,751 probes targeting 11,026 loci in two Neotropical Bignoniaceae tree species: Handroanthus serratifolius from the seasonally dry tropical forest (SDTF) and Tabebuia aurea from savannas, and compared with the population genomics of H. impetiginosus from SDTF. OutFLANK detected 29 loci in 20 genes with selection signal in H. serratifolius and no loci in T. aurea. Using BayPass, we found evidence of selection in 335 loci in 312 genes in H. serratifolius, 101 loci in 92 genes in T. aurea, and 448 loci in 416 genes in H. impetiginosus. All approaches evidenced several genes affecting plant response to environmental stress and primary metabolic processes. The three species shared no SNPs with selection signal, but we found SNPs affecting the same gene in pair of species. Handroanthus serratifolius showed differences in allele frequencies at SNPs with selection signal among ecosystems, mainly between Caatinga/Cerrado and Atlantic Forest, while H. impetiginosus had one allele fixed across all populations, and T. aurea had similar allele frequency distribution among ecosystems and polymorphism across populations. Taken together, our results indicate that natural selection related to environmental stress shaped the spatial pattern of genetic diversity in the three species. However, the three species have different geographical distribution and niches, which may affect tolerances and adaption, and natural selection may lead to different signatures due to the differences in adaptive landscapes in different niches.

RNA-Seq reveals different responses to drought in Neotropical trees from savannas and seasonally dry forests.

BACKGROUND: Water is one of the main limiting factors for plant growth and crop productivity. Plants constantly monitor water availability and can rapidly adjust their metabolism by altering gene expression. This leads to phenotypic plasticity, which aids rapid adaptation to climate changes. Here, we address phenotypic plasticity under drought stress by analyzing differentially expressed genes (DEG) in four phylogenetically related neotropical Bignoniaceae tree species: two from savanna, Handroanthus ochraceus and Tabebuia aurea, and two from seasonally dry tropical forests (SDTF), Handroanthus impetiginosus and Handroanthus serratifolius. To the best of our knowledge, this is the first report of an RNA-Seq study comparing tree species from seasonally dry tropical forest and savanna ecosystems. RESULTS: Using a completely randomized block design with 4 species × 2 treatments (drought and wet) × 3 blocks (24 plants) and an RNA-seq approach, we detected a higher number of DEGs between treatments for the SDTF species H. serratifolius (3153 up-regulated and 2821 down-regulated under drought) and H. impetiginosus (332 and 207), than for the savanna species. H. ochraceus showed the lowest number of DEGs, with only five up and nine down-regulated genes, while T. aurea exhibited 242 up- and 96 down-regulated genes. The number of shared DEGs among species was not related to habitat of origin or phylogenetic relationship, since both T. aurea and H impetiginosus shared a similar number of DEGs with H. serratifolius. All four species shared a low number of enriched gene ontology (GO) terms and, in general, exhibited different mechanisms of response to water deficit. We also found 175 down-regulated and 255 up-regulated transcription factors from several families, indicating the importance of these master regulators in drought response. CONCLUSION: Our findings show that phylogenetically related species may respond differently at gene expression level to drought stress. Savanna species seem to be less responsive to drought at the transcriptional level, likely due to morphological and anatomical adaptations to seasonal drought. The species with the largest geographic range and widest edaphic-climatic niche, H. serratifolius, was the most responsive, exhibiting the highest number of DEG and up- and down-regulated transcription factors (TF).

In situ radiation explains the frequency of dioecious palms on islands.

BACKGROUND AND AIMS: Dioecy has evolved up to 5000 times in angiosperms, despite the potentially high intrinsic costs to unisexuality. Dioecy prevents inbreeding, which is especially relevant on isolated islands when gene pools are small. Dioecy is also associated with certain dispersal traits, such as fruit size and type. However, the influence of dioecy on other life history traits and island distribution remains poorly understood. Here, we test the effect of dioecy on palm (Arecaceae) speciation rates, fruit size and frequency on islands. METHODS: We used phylogenetic comparative methods to estimate the ancestral state of the sexual system and its impact on speciation rates and fruit size. Frequency of sexual systems, effect of insularity on the probability of being dioecious, and phylogenetic clustering of island dioecious vs. mainland species were inferred. Lastly, we determined the interplay of insularity and sexual system on speciation rates. KEY RESULTS: Palms repeatedly evolved different sexual systems (dioecy, monoecy and polygamy) from a hermaphrodite origin. Differences in speciation rates and fruit size among the different sexual systems were not identified. An effect of islands on the probability of the palms being dioecious was also not found. However, we found a high frequency and phylogenetic clustering of dioecious palms on islands, which were not correlated with higher speciation rates. CONCLUSIONS: The high frequency and phylogenetic clustering may be the result of in situ radiation and suggest an 'island effect' for dioecious palms, which was not explained by differential speciation rates. This island effect also cannot be attributed to long-distance dispersal due to the lack of fruit size difference among sexual systems, and particularly because palm dispersal to islands is highly constrained by the interaction between the sizes of fruit and frugivores. Taken together, we suggest that trait flexibility in sexual system evolution and the in situ radiation of dioecious lineages are the underlying causes of the outstanding distribution of palms on islands.

Higher evolutionary rates in life-history traits in insular than in mainland palms.

Isolated islands, due to the reduced interspecific competition compared to mainland habitats, present ecological opportunities for colonizing lineages. As a consequence, island lineages may be expected to experience higher rates of trait evolution than mainland lineages. However, island effects on key life-history traits of vascular plants remain underexplored at broad spatiotemporal scales, even for emblematic island clades such as palms. Here, we used phylogenetic comparative methods to evaluate potential differences in size and macroevolutionary patterns of height and fruit diameter among mainland, continental, and volcanic island palms. Further, phylogenetic beta-diversity was used to determine if lineage turnover supported an adaptive radiation scenario on volcanic islands. Volcanic island palms were taller than their continental island and mainland counterparts, whereas continental island palms exhibited smaller fruit size. Height and fruit size of palms evolved under evolutionary constraints towards an optimal value. However, scenarios of adaptive radiation and niche conservatism were not supported for the height and fruit size of volcanic and mainland palm clades, respectively, as expected. Instead, continental island palms exhibited higher evolutionary rates for height and fruit size. Insular palm assemblages (continental and volcanic) are composed of unique lineages. Beyond representing evolutionary sources of new palm lineages, our results demonstrate that insular habitats are important in shaping palm trait diversity. Also, the higher phenotypic evolutionary rates of continental island palms suggest disparate selection pressures on this habitat type, which can be an important driver of trait diversification over time. Taken together, these results stress the importance of insular habitats for conservation of functional, phylogenetic, and taxonomic diversity of palms.

Selective Sweeps Lead to Evolutionary Success in an Amazonian Hyperdominant Palm.

Despite the global importance of tropical ecosystems, few studies have identified how natural selection has shaped their megadiversity. Here, we test for the role of adaptation in the evolutionary success of the widespread, highly abundant Neotropical palm Mauritia flexuosa. We used a genome scan framework, sampling 16,262 single-nucleotide polymorphisms (SNPs) with target sequence capture in 264 individuals from 22 populations in rainforest and savanna ecosystems. We identified outlier loci as well as signal of adaptation using Bayesian correlations of allele frequency with environmental variables and detected both selective sweeps and genetic hitchhiking events. Functional annotation of SNPs with selection footprints identified loci affecting genes related to adaptation to environmental stress, plant development, and primary metabolic processes. The strong differences in climatic and soil variables between ecosystems matched the high differentiation and low admixture in population Bayesian clustering. Further, we found only small differences in allele frequency distribution in loci putatively under selection among widespread populations from different ecosystems, with fixation of a single allele in most populations. Taken together, our results indicate that adaptive selective sweeps related to environmental stress shaped the spatial pattern of genetic diversity in M. flexuosa, leading to high similarity in allele frequency among populations from different ecosystems.

Megafauna Seed Dispersal in the Neotropics: A Meta-Analysis Shows No Genetic Signal of Loss of Long-Distance Seed Dispersal.

Restricted gene flow may lead to the loss of genetic diversity and higher genetic differentiation among populations, but the genetic consequences of megafauna extinction for plant populations still remain to be assessed. We performed a phylogenetic-independent meta-analysis across 102 Neotropical plants to test the hypothesis that plant species with megafaunal seed dispersal syndrome have a lower genetic diversity and a higher genetic differentiation than those without it. We classified as megafauna-dependent plant species those that potentially relied only on megafauna to seed dispersal, and as megafauna-independent those that relied on megafauna and other seed dispersers. Our data comprised 98 studies using microsatellite markers. We found no statistical difference in genetic diversity and differentiation between plants with megafauna and non-megafauna seed dispersal syndrome, although the statistical power to detect differences in genetic differentiation was low. Moreover, we found no statistical difference between megafauna-dependent and megafauna-independent plant species. We then used generalized linear mixed models and phylogenetic generalized least square models to investigate the effects of megafaunal seed dispersal syndromes and reproductive traits on variation in genetic diversity and genetic differentiation. We found no effect of megafaunal syndrome, rather, reproductive traits, such as pollination mode, mating, and breeding systems, showed significant effects. Our findings show that the genetic studies of Neotropical plants performed so far show no difference in genetic diversity and differentiation in plants with megafaunal compared to those with non-megafaunal seed dispersal syndromes. Our results also provide evidence pointing out that plant species with megafaunal seed dispersal syndromes may have used different strategies to counterbalance the extinction of their mutualistic megafauna dispersers, such as the dispersal by extant mammals that may promote long-distance seed dispersal. Our results also reinforce the importance of pollination to long-distance gene flow in Neotropical plants.

Adjacency and Area Explain Species Bioregional Shifts in Neotropical Palms.

Environmental and geographical variables are known drivers of community assembly, however their influence on phylogenetic structure and phylogenetic beta diversity of lineages within different bioregions is not well-understood. Using Neotropical palms as a model, we investigate how environmental and geographical variables affect the assembly of lineages into bioregions across an evolutionary time scale. We also determine lineage shifts between tropical (TRF) and non-tropical (non-TRF) forests. Our results identify that distance and area explain phylogenetic dissimilarity among bioregions. Lineages in smaller bioregions are a subset of larger bioregions and contribute significantly to the nestedness component of phylogenetic dissimilarity, here interpreted as evidence for a bioregional shift. We found a significant tendency of habitat shifts occurring preferentially between TRF and non-TRF bioregions (31 shifts) than from non-TRF to TRF (24) or from TRF to TRF (11) and non-TRF to non-TRF (9). Our results also present cases where low dissimilarity is found between TRF and non-TRF bioregions. Most bioregions showed phylogenetic clustering and larger bioregions tended to be more clustered than smaller ones, with a higher species turnover component of phylogenetic dissimilarity. However, phylogenetic structure did not differ between TRF and non-TRF bioregions and diversification rates were higher in only two lineages, Attaleinae and Bactridinae, which are widespread and overabundant in both TRF and non-TRF bioregions. Area and distance significantly affected Neotropical palm community assembly and contributed more than environmental variables. Despite palms being emblematic humid forest elements, we found multiple shifts from humid to dry bioregions, showing that palms are also important components of these environments.

A genome-wide scan shows evidence for local adaptation in a widespread keystone Neotropical forest tree.

The role of natural selection in shaping patterns of diversity is still poorly understood in the Neotropics. We carried out the first genome-wide population genomics study in a Neotropical tree, Handroanthus impetiginosus (Bignoniaceae), sampling 75,838 SNPs by sequence capture in 128 individuals across 13 populations. We found evidences for local adaptation using Bayesian correlations of allele frequency and environmental variables (32 loci in 27 genes) complemented by an analysis of selective sweeps and genetic hitchhiking events using SweepFinder2 (81 loci in 47 genes). Fifteen genes were identified by both approaches. By accounting for population genetic structure, we also found 14 loci with selection signal in a STRUCTURE-defined lineage comprising individuals from five populations, using Outflank. All approaches pinpointed highly diverse and structurally conserved genes affecting plant development and primary metabolic processes. Spatial interpolation forecasted differences in the expected allele frequencies at loci under selection over time, suggesting that H. impetiginosus may track its habitat during climate changes. However, local adaptation through natural selection may also take place, allowing species persistence due to niche evolution. A high genetic differentiation was seen among the H. impetiginosus populations, which, together with the limited power of the experiment, constrains the improved detection of other types of soft selective forces, such as background, balanced, and purifying selection. Small differences in allele frequency distribution among widespread populations and the low number of loci with detectable adaptive sweeps advocate for a polygenic model of adaptation involving a potentially large number of small genome-wide effects.

Design and evaluation of a sequence capture system for genome-wide SNP genotyping in highly heterozygous plant genomes: a case study with a keystone Neotropical hardwood tree genome.

Targeted sequence capture coupled to high-throughput sequencing has become a powerful method for the study of genome-wide sequence variation. Following our recent development of a genome assembly for the Pink Ipê tree (Handroanthus impetiginosus), a widely distributed Neotropical timber species, we now report the development of a set of 24,751 capture probes for single-nucleotide polymorphisms (SNPs) characterization and genotyping across 18,216 distinct loci, sampling more than 10 Mbp of the species genome. This system identifies nearly 200,000 SNPs located inside or in close proximity to almost 14,000 annotated protein-coding genes, generating quality genotypic data in populations spanning wide geographic distances across the species native range. To provide recommendations for future developments of similar systems for highly heterozygous plant genomes we investigated issues such as probe design, sequencing coverage and bioinformatics, including the evaluation of the capture efficiency and a reassessment of the technical reproducibility of the assay for SNPs recall and genotyping precision. Our results highlight the value of a detailed probe screening on a preliminary genome assembly to produce reliable data for downstream genetic studies. This work should inspire and assist the development of similar genomic resources for other orphan crops and forest trees with highly heterozygous genomes.

Unravelling the genetic differentiation among varieties of the Neotropical savanna tree Hancornia speciosa Gomes.

Background and Aims: Spatial distribution of species genetic diversity is often driven by geographical distance (isolation by distance) or environmental conditions (isolation by environment), especially under climate change scenarios such as Quaternary glaciations. Here, we used coalescent analyses coupled with ecological niche modelling (ENM), spatially explicit quantile regression analyses and the multiple matrix regression with randomization (MMRR) approach to unravel the patterns of genetic differentiation in the widely distributed Neotropical savanna tree, Hancornia speciosa (Apocynaceae). Due to its high morphological differentiation, the species was originally classified into six botanical varieties by Monachino, and has recently been recognized as only two varieties by Flora do Brasil 2020. Thus, H. speciosa is a good biological model for learning about evolution of phenotypic plasticity under genetic and ecological effects, and predicting their responses to changing environmental conditions. Methods: We sampled 28 populations (777 individuals) of Monachino's four varieties of H. speciosa and used seven microsatellite loci to genotype them. Key Results: Bayesian clustering showed five distinct genetic groups (K = 5) with high admixture among Monachino's varieties, mainly among populations in the central area of the species geographical range. Genetic differentiation among Monachino's varieties was lower than the genetic differentiation among populations within varieties, with higher within-population inbreeding. A high historical connectivity among populations of the central Cerrado shown by coalescent analyses may explain the high admixture among varieties. In addition, areas of higher climatic suitability also presented higher genetic diversity in such a way that the wide historical refugium across central Brazil might have promoted the long-term connectivity among populations. Yet, FST was significantly related to geographic distances, but not to environmental distances, and coalescent analyses and ENM predicted a demographical scenario of quasi-stability through time. Conclusions: Our findings show that demographical history and isolation by distance, but not isolation by environment, drove genetic differentiation of populations. Finally, the genetic clusters do not support the two recently recognized botanical varieties of H. speciosa, but partially support Monachino's classification at least for the four sampled varieties, similar to morphological variation.

The road to evolutionary success: insights from the demographic history of an Amazonian palm.

Evolutionary success, as demonstrated by high abundance and a wide geographical range, is related to genetic variation and historical demography. Here we assess how climatic change during the Quaternary influenced the demography and distribution of the Neotropical swamp palm Mauritia flexuosa. Using microsatellite loci and coalescent analyses we examined how demographical dynamics affected genetic diversity, effective population size and connectivity through time and space. Mauritia flexuosa presents significant genetic differentiation between the Amazonian and Cerrado biomes and among different river basins. Amazonian lineages are ancient compared to lineages from the Cerrado, a pattern corroborated using the fossil pollen record, where the species was absent from the Cerrado during the cold and dry periods of the last glacial cycles, then returned during the wet, interglacial phases. Coalescent simulations show that the pattern of observed genetic diversity for M. flexuosa is most likely due to a range retraction during the Last Glacial Maximum, leading to multiple refugia and resulting in high differentiation between Amazonian and Cerrado biomes. Isolation-by-distance and by-environment also shaped the distribution and evolutionary success of M. flexuosa. Our study provides new insights into the historical factors that affected geographical distribution and structure genetic diversity, contributing to long-term evolutionary success.

Genome assembly of the Pink Ipê (Handroanthus impetiginosus, Bignoniaceae), a highly valued, ecologically keystone Neotropical timber forest tree.

Background: Handroanthus impetiginosus (Mart. ex DC.) Mattos is a keystone Neotropical hardwood tree widely distributed in seasonally dry tropical forests of South and Mesoamerica. Regarded as the "new mahogany," it is the second most expensive timber, the most logged species in Brazil, and currently under significant illegal trading pressure. The plant produces large amounts of quinoids, specialized metabolites with documented antitumorous and antibiotic effects. The development of genomic resources is needed to better understand and conserve the diversity of the species, to empower forensic identification of the origin of timber, and to identify genes for important metabolic compounds. Findings: The genome assembly covers 503.7 Mb (N50 = 81 316 bp), 90.4% of the 557-Mbp genome, with 13 206 scaffolds. A repeat database with 1508 sequences was developed, allowing masking of ∼31% of the assembly. Depth of coverage indicated that consensus determination adequately removed haplotypes assembled separately due to the extensive heterozygosity of the species. Automatic gene prediction provided 31 688 structures and 35 479 messenger RNA transcripts, while external evidence supported a well-curated set of 28 603 high-confidence models (90% of total). Finally, we used the genomic sequence and the comprehensive gene content annotation to identify genes related to the production of specialized metabolites. Conclusions: This genome assembly is the first well-curated resource for a Neotropical forest tree and the first one for a member of the Bignoniaceae family, opening exceptional opportunities to empower molecular, phytochemical, and breeding studies. This work should inspire the development of similar genomic resources for the largely neglected forest trees of the mega-diverse tropical biomes.

Correction: Pollination Mode and Mating System Explain Patterns in Genetic Differentiation in Neotropical Plants.

[This corrects the article DOI: 10.1371/journal.pone.0158660.].

Climatic changes can drive the loss of genetic diversity in a Neotropical savanna tree species.

The high rates of future climatic changes, compared with the rates reported for past changes, may hamper species adaptation to new climates or the tracking of suitable conditions, resulting in significant loss of genetic diversity. Trees are dominant species in many biomes and because they are long-lived, they may not be able to cope with ongoing climatic changes. Here, we coupled ecological niche modelling (ENM) and genetic simulations to forecast the effects of climatic changes on the genetic diversity and the structure of genetic clusters. Genetic simulations were conditioned to climatic variables and restricted to plant dispersal and establishment. We used a Neotropical savanna tree as species model that shows a preference for hot and drier climates, but with low temperature seasonality. The ENM predicts a decreasing range size along the more severe future climatic scenario. Additionally, genetic diversity and allelic richness also decrease with range retraction and climatic genetic clusters are lost for both future scenarios, which will lead genetic variability to homogenize throughout the landscape. Besides, climatic genetic clusters will spatially reconfigure on the landscape following displacements of climatic conditions. Our findings indicate that climate change effects will challenge population adaptation to new environmental conditions because of the displacement of genetic ancestry clusters from their optimal conditions.

Population structure and genetic diversity of the giant anteater (Myrmecophaga tridactyla: Myrmecophagidae, Pilosa) in Brazil

Abstract The giant anteater (Myrmecophaga tridactyla, Pilosa, Linnaeus 1758) belongs to the mammalian order Pilosa and presents a large distribution along South America, occupying a great variety of habitats. It is listed in the IUCN Red List of threatened species as Vulnerable. Despite threatened, there is a lack of studies regarding its genetic variability. The aim of this study was to examine the genetic diversity and patterns of genetic structure within remaining populations. We analyzed 77 individuals from seven different populations distributed in four biomes across Brazil: Cerrado, Pantanal, Atlantic Forest and Amazon Forest. We sequenced two mitochondrial markers (control region and Cyt-b) and two nuclear markers (AMELY and RAG2). We found high genetic diversity within subpopulations from National Parks of Serra da Canastra and Emas, both within the Cerrado biome, with signs of population expansion. Besides, we found a notable population structure between populations from the Cerrado/Pantanal and Amazon Forest biomes. This data is a major contribution to the knowledge of the evolutionary history of the species and to future management actions concerning its conservation.

Population structure and genetic diversity of the giant anteater (Myrmecophaga tridactyla: Myrmecophagidae, Pilosa) in Brazil.

The giant anteater (Myrmecophaga tridactyla, Pilosa, Linnaeus 1758) belongs to the mammalian order Pilosa and presents a large distribution along South America, occupying a great variety of habitats. It is listed in the IUCN Red List of threatened species as Vulnerable. Despite threatened, there is a lack of studies regarding its genetic variability. The aim of this study was to examine the genetic diversity and patterns of genetic structure within remaining populations. We analyzed 77 individuals from seven different populations distributed in four biomes across Brazil: Cerrado, Pantanal, Atlantic Forest and Amazon Forest. We sequenced two mitochondrial markers (control region and Cyt-b) and two nuclear markers (AMELY and RAG2). We found high genetic diversity within subpopulations from National Parks of Serra da Canastra and Emas, both within the Cerrado biome, with signs of population expansion. Besides, we found a notable population structure between populations from the Cerrado/Pantanal and Amazon Forest biomes. This data is a major contribution to the knowledge of the evolutionary history of the species and to future management actions concerning its conservation.

Demographic stability and high historical connectivity explain the diversity of a savanna tree species in the Quaternary.

Background and Aims: Cyclic glaciations were frequent throughout the Quaternary and this affected species distribution and population differentiation worldwide. The present study reconstructed the demographic history and dispersal routes of Eugenia dysenterica lineages and investigated the effects of Quaternary climate change on its spatial pattern of genetic diversity. Methods: A total of 333 individuals were sampled from 23 populations and analysed by sequencing four regions of the chloroplast DNA and the internal transcribed spacer of the nuclear DNA. The analyses were performed using a multi-model inference approach based on ecological niche modelling and statistical phylogeography. Key Results: Coalescent simulation showed that population stability through time is the most likely scenario. The palaeodistribution dynamics predicted by the ecological niche models revealed that the species was potentially distributed across a large area, extending over Central-Western Brazil through the last glaciation. The lineages of E. dysenterica dispersed from Central Brazil towards populations at the northern, western and south-eastern regions. A historical refugium through time may have favoured lineage dispersal and the maintenance of genetic diversity. Conclusions: The results suggest that the central region of the Cerrado biome is probably the centre of distribution of E. dysenterica and that the spatial pattern of its genetic diversity may be the outcome of population stability throughout the Quaternary. The lower genetic diversity in populations in the south-eastern Cerrado biome is probably due to local climatic instability during the Quaternary.

Demographical history and palaeodistribution modelling show range shift towards Amazon Basin for a Neotropical tree species in the LGM.

BACKGROUND: We studied the phylogeography and demographical history of Tabebuia serratifolia (Bignoniaceae) to understand the disjunct geographical distribution of South American seasonally dry tropical forests (SDTFs). We specifically tested if the multiple and isolated patches of SDTFs are current climatic relicts of a widespread and continuously distributed dry forest during the last glacial maximum (LGM), the so called South American dry forest refugia hypothesis, using ecological niche modelling (ENM) and statistical phylogeography. We sampled 235 individuals of T. serratifolia in 17 populations in Brazil and analysed the polymorphisms at three intergenic chloroplast regions and ITS nuclear ribosomal DNA. RESULTS: Coalescent analyses showed a demographical expansion at the last c. 130 ka (thousand years before present). Simulations and ENM also showed that the current spatial pattern of genetic diversity is most likely due to a scenario of range expansion and range shift towards the Amazon Basin during the colder and arid climatic conditions associated with the LGM, matching the expected for the South American dry forest refugia hypothesis, although contrasting to the Pleistocene Arc hypothesis. Populations in more stable areas or with higher suitability through time showed higher genetic diversity. Postglacial range shift towards the Southeast and Atlantic coast may have led to spatial genome assortment due to leading edge colonization as the species tracks suitable environments, leading to lower genetic diversity in populations at higher distance from the distribution centroid at 21 ka. CONCLUSION: Haplotype sharing or common ancestry among populations from Caatinga in Northeast Brazil, Atlantic Forest in Southeast and Cerrado biome and ENM evince the past connection among these biomes.